Percolation transition and dissipation in quantum Ising magnets 
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We study the effects of dissipation on a randomly diluted transverse-field Ising magnet close to 
the percolation threshold. For weak transverse fields, a novel percolation quantum phase transition 
separates a superparamagnetic cluster phase from an inhomogeneously ordered ferromagnetic phase. 
The properties of this transition are dominated by large frozen and slowly fluctuating percolation 
clusters. This leads to a discontinuous magnetization- field curve and exotic hysteresis phenomena 
as well as highly singular behavior of magnetic susceptibility and specific heat. We compare our 
results to the smeared transition in generic dissipative random quantum Ising magnets. We also 
discuss the relation to metallic quantum magnets and other experimental realizations. 



In diluted quantum magnets, the combination of ge- 
ometric and quantum fluctuations can lead to uncon- 
ventional low-temperature properties such as singular 
thermodynamic quantities in quantum Griffiths phases 
as well as exotic scaling at the quantum phase transi- 
tions (QPTs) between magnetic and nonmagnetic ground 
states 1-4 (for a recent review see, e.g., Ref. 5). In many 



real systems, the magnetic degrees of freedom are coupled 
to an environment of "heat bath" modes. The dissipa- 
tion caused by the bath is known to qualitatively change 
the properties even of a single quantum spin. 6-9 . 

In this Letter, we show that dissipation dramatically 
changes phases and phase transitions of a randomly di- 
luted quantum Ising magnet, as is illustrated by Fig. 1. 
The behavior without dissipation is well understood: 10-13 
Magnetic long-range order exists for dilutions p up to the 
lattice percolation threshold p c as long as the transverse 
field h x is below a critical field h c (p). For p > p c , long- 
range order is impossible because the system is decom- 
posed into noninteracting finite-size clusters. There are 
two QPTs, separated by a multicritical point at (p c , h*). 
The transition across h c (j>) for < p < p c (transition 
(i) in Fig. la) falls into the generic random quantum 
Ising universality class and is characterized by an infinite- 
randomness fixed point. 1,2 The transition across p c at 
h x < h* (transition (ii) in Fig. la) is an unusual kind of 
percolation QPT. 13 The paramagnetic phase consists of 
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FIG. 1: Schematic ground state phase diagrams of diluted 
quantum Ising magnets without (a) and with (b) dissipation. 
CSPM is the classical superparamagnetic phase. The dashed 
line in (b) marks the crossover between homogeneous and 
inhomogeneous order in the smeared transition scenario. 



two regions: a conventional gapped region for h x > h c (0) 
and a gapless quantum Griffiths region for h x < /i c (0) 
where rare regions lead to power-law thermodynamic sin- 
gularities. 

In the presence of weak ohmic dissipation, the phase 
diagram changes qualitatively since the quantum dynam- 
ics of sufficiently large clusters freezes. 14-16 This leads 
to the appearance of an unusual classical superparamag- 
netic cluster phase. For p > p c , it is stabilized down to 
zero temperature because different clusters are strictly 
decoupled. In contrast, for p < p c weakly coupled frozen 
rare regions align. Magnetic long-range order thus ex- 
tends to the clean critical field h CjCt , and the field-driven 
phase transition (i) is smeared. 17-19 

In this Letter, we focus on the novel percolation QPT 
at small fields h x (transition (ii) in Fig. lb). Our results 
can be summarized as follows: The total magnetization 
consists of a coherent part TOqo and an incoherent part 
m s t . The critical behavior of the coherent part is given by 
classical percolation theory ~ \p — p c \^ a with j3 c the 
classical percolation exponent. Somewhat unexpectedly, 
moo + TUst (which is measured in an infinitesimal order- 
ing field) is noncritical at p c . The interplay between the 
two parts leads to unusual magnetization-field hysteresis. 
At low temperatures T, the susceptibility has a Curie 
contribution from the frozen clusters that diverges as 
X ~ \p~Pc\~ lc /T with the classical percolation exponent 
7 C . Fluctuating clusters provide a subleading l/(Tln 2 T) 
contribution. In contrast, the low-temperature specific 
heat is dominated by fluctuating clusters, resulting in 
a logarithmic temperature dependence, C ~ ln -2 (l/T). 
The logarithmic terms in \ an d C exist on both sides of 
p c and their prefactors are noncritical at this point . 

Our starting point is a d-dimensional (d > 2) site- 
diluted transverse-field Ising model 10-13 
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a prototypical disordered quantum magnet. The Pauli 
matrices of and af represent the spin components at site 
i, the exchange interaction J couples nearest neighbor 
sites, and the transverse field h x controls the quantum 
fluctuations. Dilution is introduced via random variables 
which can take the values and 1 with probabilities 
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p and 1 — p, respectively. We now couple each spin to a 
local ohmic bath of harmonic oscillators, 19,20 



1 — sa < A s /uj c , eq. (3) has to be replaced by ; 



2H 



, (2) 



where a^n and a ■ n arc the annihilation and creation op- 
erators of the n-th oscillator coupled to spin i; v>i tTl is its 
frequency, and Xi iTl is the coupling constant. All baths 
have the same spectral function £(ui) — tt J2 n Xf n 5(ui — 

v i,n)l v i,n = 27rau;e _a with a the dimensionlcss dis- 
sipation strength and co c the cutoff energy. Such local 
ohmic dissipation can be realized in good approximation 
in magnetic nanoparticles in an insulating host. 21 It also 
applies to resistively shunted SQUIDs 6 and qualitatively 
(possibly with a different spectral density £(ui)) to molec- 
ular magnets weakly coupled to nuclear spins 22 (see, e.g., 
the spin- 1/2 molecular complex V15, Ref. 23). 

Let us start by considering a single percolation clus- 
ter of s occupied sites. Without dissipation, and for 
small transverse fields, h x <C h c . a ~ J, its lowest two 
energy levels correspond to the states of a single effec- 
tive Ising spin with a moment proportional to s in an 
effective transverse field (tunneling matrix element) A s . 
For large s, A s can be estimated in s-th order of pertur- 
bation theory to be h x e~ Bs , with B a constant of order 
ln(J//i x ). 13 All other levels are separated by energies of 
order J and can thus be neglected for the low-energy 
physics. This means the cluster tunnels coherently be- 
tween up and down with a tunneling frequency A s . 

The effects of the heat bath on the cluster can be 
worked out in 2nd order perturbation theory in h x . We 
find that the effective spin of the cluster feels a sin- 
gle ohmic bath with dissipation strength Asa because 
the number of oscillators coupling to the cluster is pro- 
portional to s. 24 Here, A is an /i^-dependent constant 
which we suppress as it can be absorbed in the cluster 
moment. Analogous results have been obtained for the 
dissipative random quantum Ising chain 19 and itinerant 
Ising magnets. 15 Each percolation cluster thus behaves 
like an ohmic spin-boson model: 6,7 For strong dissipa- 
tion, sa > a c w l 37 , the cluster is in the localized phase, 
i.e., the renormalizcd tunneling matrix element A^ s van- 
ishes, and the magnetization is frozen in one of the di- 
rections. Clusters with sa < 1 are in the delocalized 
phase, i.e., they still tunnel, but with a greatly reduced 
frequency. For sa <C 1 it can be estimated by adiabatic 
renormalization 6 , 

A R , S ~ A s (A s K) as/(1 - Qs) - h x exp [-bs/(l - sa)} , 

(3) 

with the constant b = B + a\n(u c /h x ). The QPT at 
sa = 1 is of Kostcrlitz-Thoulcss 25 type, analogous to 
the transition in the classical 1/r 2 Ising chain. 26,27 Here, 
Ar. s plays the role of the classical correlation length. 
Therefore, the functional form (3) remains valid as long 
as 1 — sa > A s /uo c . However, very close to the transition, 



h x exp (— 6's/Vl — sa) 



(4) 



These results allows us to determine the phase dia- 
gram, Fig. lb, for fixed small dissipation strength a. For 
h x > h CM , our system is a conventional gapped quantum 
paramagnet for all p because not even the clean bulk 
system orders. For h x < h c _ a and p > p c , the system 
is decomposed into finite-size percolation clusters. The 
largest of these clusters are frozen and behave like classi- 
cal Ising magnets while the smaller ones fluctuate. The 
system is thus in a classical superparamagnetic cluster 
phase. For p < p c , ferromagnetic long-range order exists 
on the infinite percolation cluster. The ferromagnetic 
phase extends up to the clean critical field h c , a for all 
p < p c . This is a manifestation of the smeared transi- 
tion scenario: 17 The infinite percolation cluster contains, 
with small probability, large impurity-free regions that 
freeze close to h c a . Due to their coupling to the rest of 
the infinite cluster they induce a highly inhomogeneous 
but coherent magnetization already for transverse fields 
close to h Ct<x . For smaller transverse fields there will be 
a crossover to homogeneous order on the infinite cluster 
(see the dashed line in Fig. lb). 

To study the percolation QPT, we need the distribu- 
tion of cluster sizes as a function of dilution p. According 
to percolation theory, 29 the number n s of occupied clus- 
ters of size s per lattice site obeys the scaling form 



(5) 



Here t = p — p c , and r c and a c are classical percolation 
exponents 38 . The scaling function f(x) is analytic for 
small x and has a single maximum at some x n 
For large it drops off rapidly 



exp(-(c ia ;) 1 /^) 



(for x > 0), 



exp (-(c2x)( 1 - 1 / d V a °) (for x < 0), 



> 0. 



(6) 



where c\ and C2 are constants of order unity. In addi- 
tion, for p < p c , there is an infinite cluster containing a 
fraction ~ \p — p c \^ c of all sites. All classical perco- 
lation exponents are determined by r c and a c including 
the correlation length exponent v c = (r c — l)/(da c ), the 
order parameter exponent f3 c — (r c — 2)/<r c , and the sus- 
ceptibility exponent j c = (3 — t c )/<t c . 

The low-energy density of states (DOS) Pdy(e) of the 
dynamic clusters (sa < 1) can be determined by combin- 
ing the single-cluster results (3) and (4) with the cluster 
size distribution (5) via /fdy(e) = X) s <s n s5(e — A^ s ). 
Inserting (3) for Ar. s yields 



Pdy(e) = n s{e) b/ [e (b + a ln(h x /e)) 2 



(7) 



where 



is the density of clusters of size s(e) = 



\n{h x /e)/[b + aln(h x /e)]. Eq. (7) contains the full 
crossover between quantum Griffiths behavior at higher 
energies and damping-dominated behavior at low ener- 
gies. For e > e cross = h x e~ b / a , the DOS simplifies to a 
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FIG. 2: Zero-temperature magnetization as a function of di- 
lution p for different ordering fields H z . The solid line is the 
coherent magnetization at H z = (infinite cluster only). The 
dashed line is for an infinitesimal field. The insets show the 
singular field dependence of m below and above p c . 

power-law form, Pdy(e) ~ e c / b_1 , characteristic of quan- 
tum Griffiths behavior. 3 In contrast, for e < e cross , the 
energy dependence is even more singular, 

p iY {e)~n s J[e\^{h x /e)\ , (8) 

with (f> — 2. At very low energies, in the asymptotic re- 
gion of the Kosterlitz-Thouless transition, A^ s is given 
by (4) which leads to cf> — 3. The functional form (8) of 
the low-energy DOS is the same on both sides of the per- 
colation threshold since it is caused by clusters of finite 
size s w s c = 1/a. Moreover, the prefactor is noncrit- 
ical at p c because the cluster size distribution n s is an 
analytic function of p for any finite s. 

We now discuss the physics at the percolation tran- 
sition in small transverse fields (h x <C h CtCt ), starting 
with the total magnetization to. We have to distinguish 
the contributions md y from dynamic clusters, m s t from 
frozen finite-size clusters, and TOoo from the infinite per- 
colation cluster, if any. For zero ordering field H z in z- 
direction, md y vanishes, because the dynamic clusters are 
in symmetric superpositions of the up and down states. 
The frozen finite-size clusters individually have a nonzero 
magnetization, but it sums up to m st — because they 
do not align coherently for H z = 0. The only coherent 
contribution to the total magnetization is m^. Since the 
infinite cluster is long-range ordered for small transverse 
field h Xl its magnetization is proportional to the number 
Poo of sites in the infinite cluster, giving 

m.m.~P. W ~{l'-*l* fZltli W 

(solid line in Fig. 2) . The magnetization critical exponent 
is given by its classical lattice percolation value (3 C 

The zero-temperature response of the magnetization to 
a small ordering field H z is highly singular. Frozen finite- 
size clusters align parallel to the field already for infinites- 
imal H Z) leading to a jump in m(H z ) at H z = 0. The 
magnitude of the jump is given by m st ~ X) s >s n s- At 



the percolation threshold, m st is approximately given by 
(1 — p c )s 2_Tc , and it is exponentially small for both p — > 
and p — > 1. Somewhat surprisingly, the total magnetiza- 
tion in an infinitesimal field, given by TOqq + m st , is ana- 
lytic (i.e., noncritical) at p = p c (see dashed line in Fig. 
2). This follows from the fact that only clusters with sizes 
below s c are not polarized. Since the cluster size distribu- 
tion is analytic for any finite s, their total size is analytic 
at p c , and consequently TOoo +TO s t must be analytic, too. 
To estimate the contribution mdy of the dynamic clusters 
to the magnetization, we integrate the magnetization of 
a single cluster, H z s 2 (e)/(H 2 s 2 (e) + e 2 ) 1 / 2 over the DOS 
given in (7). For small fields H z < e cross /s c , the main 
contribution comes from clusters with lowest energies and 
sizes close to s c , resulting in a highly singular dependence 
of mdy on the ordering field, TOd y ~ n Sc / \d^~ 1 {h x / H z ) 
for H z — > 0. For higher fields, we find a crossover to 

c/b 

quantum Griffiths behavior of the form TOdy ~ H z . 

As a consequence, our system displays unusual hys- 
teresis (insets of Fig. 2). For p < p c , the long-range 
ordered infinite cluster contributes a conventional single- 
domain hysteresis loop. Finite-size frozen clusters do not 
show hysteresis, for arbitrarily small H z , they relax with 
a nonzero rate. 6 Thus, they contribute jumps in m at 
H z — 0. The dynamic clusters add logarithmic singulari- 
ties at H z = 0. For p > p c , there is no infinite cluster, but 
the jump at H z = as well as the logarithmic singularity 
survive. The above picture is valid for adiabatically slow 
sweeps of the ordering field H z . If H z is changed at a 
finite rate, the largest clusters will fall out of equilibrium 
which leads to very interesting time-dependent hysteresis 
phenomena that will be reported elsewhere. 24 

The low-tempcrature susceptibility is dominated by 
the contribution Xst of the static clusters, each of which 
has a Curie susceptibility of the form s 2 /T. Summing 
over all static clusters gives 

Xst ~ n * s2 l T ~ \P ~ Pc\~ 1C / T ( 10 ) 

s>s c 

close to the percolation threshold. For p — * and 
p — > 1, the prefactor of the Curie term vanishes ex- 
ponentially. To determine the subleading contribution 
Xdy of the dynamic clusters, we integrate the suscep- 
tibility [s 2 (e) tanh(e/T)]/e of a single cluster over the 
DOS (7). For T < e cross , %dy is dominated by the 
clusters with lowest energies and sizes close to s c , giv- 
ing Xdy ~ n Sa /[Thi^~' 1 (h x /T)]. For temperatures above 
Across, damping becomes unimportant and Xdy takes the 
quantum Griffiths form Xdy ~ j^c/b-i^ 

Finally, we consider the heat capacity which is domi- 
nated by finite size clusters, the infinite cluster does not 
contribute for T <C J. The heat capacity of a single 
dynamic cluster with gap e is not equal to that of a two- 
level system, but given by C w 7ras(e)T/(3e) for T <C e 
and C ~ (e/T) 2 - 2s W Q for T > e. 30 Integrating over 
the DOS (8) gives C dy - n s J In* -1 (h x /T). The frozen 
clusters contribute a term of the same functional form, 
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C st ~ n Be /\n! l '- 1 (h x /T). For T > e cross , the logarith- 
mic temperature dependence is replaced by the quantum 
Griffiths behavior C ~ T c / b . 

To summarize, we have have shown that the diluted 
quantum Ising model with ohmic dissipation undergoes 
a novel percolation QPT. Observables dominated by the 
frozen clusters, such as total magnetization and suscep- 
tibility, show classical percolation critical behavior. For 
these quantities our analysis is asymptotically exact. Ob- 
servables dominated by the fluctuating clusters, such as 
the heat capacity, are noncritical across the percolation 
threshold but display singular temperature dependencies 
(that can be expressed in terms of the low-energy prop- 
erties of the ohmic spin-boson model) . In the remaining 
paragraphs we put our results into broader perspective. 

In contrast to the generic transition for p < p c (transi- 
tion (i) in Fig. lb), our percolation transition (transition 
(ii) in Fig. lb) is not smeared by the mechanism of Ref. 
17 because different percolation clusters are not coupled. 
Instead, the frozen clusters make an incoherent contribu- 
tion to the magnetization. Deviations from a pure perco- 
lation scenario can change this. If tails in the interaction 
couple (even very weakly) different frozen clusters, they 
align coherently, and the smeared transition scenario is 
restored. However, for a rapidly decaying interaction, the 
smearing of the percolation transition becomes important 
only at very low energies. 



In addition to the experimental examples mentioned 
after (2), our dissipative quantum Ising model is related 
to itinerant Ising magnets where the damping is due to 
the electrons. Integrating out the bath modes in (2) leads 
to a model in the same universality class as Hertz' 31 field 
theory for the itinerant antiferromagnetic transition. 20 ' 32 
Thus, local ohmic dissipation correctly captures the lead- 
ing effects of the damping of the magnetic modes in 
itinerant magnets. It should be emphasized, however, 
that itinerant magnets contain extra complications such 
as long-range RKKY interactions, making a pure perco- 
lation scenario very unlikely. In Ref. 16, a theory simi- 
lar to our description of the classical superparamagnetic 
phase was suggested as an explanation of the unusual 
thermodynamics 33 of heavy fermion compounds. Such a 
theory can be expected to hold above the smeared phase 
transition. 17 Moreover, the importance of the dissipation 
in these systems is currently controversial. 34-36 

Other possible applications include two-level atoms in 
diluted optical lattices coupled to an electromagnetic 
field, random arrays of tunneling defects in solids or, in 
the future, many coupled qubits in a noisy environment. 
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